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Abstract. We consider the set Q2 of double zeros in (0, 1) for power series with 
£ — ' coefficients in { — 1, 0, I}. We prove that Q.2 is disconnected, and estimate minf22 

with high accuracy. We also show that [2 1//2 — 77, I) C SI2 for some small, but 
explicit 77 > (this was only known for 77 = 0). These results have applications 
' in the study of infinite Bernoulli convolutions and connectedness properties of 

self-affine fractals. 

•5 
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1. Introduction 



Let 



71=1 



in! B={l + Y^a n x n : on €{-1,0,1}}. (1.1) 

IT) 

■^j- ■ We investigate the set 

o ' 

g; n 2 = {xe(0,l) : 3/€B, f(x) = f'(x) = 0}. (1.2) 

' Thus, 0,2 is the set of zeros of power series with coefficients in { — 1,0, 1} of order 

greater or equal to two (we call them "double zeros" for short). Since B is a normal 
family, O2 is relatively closed in (0, 1), so there exists 02 = minf^- 

We prove that the set ^2 is disconnected and show that 02 ~ .6684756 (see 
Theorem 12. 5|) , In fact, we found 58 distinct components in O2, and we conjecture 
that there are infinitely many components. Numerical evidence indicates that the 
structure of Q2 is very complicated. 

Let 0L2 := sup((0, 1) \ Q 2 )- It is known that [2 _1 / 2 ,1) C 2 (see Lemma l2~2f). 
hence a 2 < 2" 1 / 2 . We show that 5 2 < 2" 1 / 2 - 4 x 10~ 6 (see Theorem l2~T0|) . 
Numerical evidence suggests that S?2 ~ 0.67, but this is harder to prove rigorously. 

The study of 0,2 is motivated by the work on infinite Bernoulli convolutions 
[TH ITU1 ITTl |5] and on some fractal sets ^5] , where a key step is checking a certain 
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transversality condition. This condition holds precisely on (0, 1) \ Q2 • Those papers 
used the estimate «2 > 0.649... obtained in ^1] by computing the smallest double 
zero of a larger class of power series B := {1 + Y^=\ a nX n ■ a n G [—1, 1]}. Theo- 
rem extends considerably (by more than 12%) the set of parameters where the 
results of [HI El d apply. 

Another motivation comes from the study of connectedness loci for certain fami- 
lies of self-similar and self-affme fractals in the plane. One of them is the "Mandel- 
brot set for pairs of linear maps" A4, studied in [31 [21 El El EEHI ED EEH] • We introduce 
two other connectedness loci, N and 0; the latter one, associated with a linear 
map having a 2 x 2 Jordan block, coincides with VL2 U (— O2). Theorem 12.51 yields 
new information about "spikes," or "antennas" — peculiar features of A4 and M. 
Our second main result is Theorem 12 . 1UI in which we obtain explicit neighborhoods 
of (previously unknown) interior points of all three connectedness loci. 

Let us make a few comments about the proofs. We use a C++ program, based 
on a modification of Bandt's algorithm from pQ, with rigorous estimates, to rule out 
double zeros in specific intervals. This program also indicates if there is a possible 
root in the interval, and provides a polynomial which is the initial part of a power 
series in B with a double zero in the interval. Once such a polynomial is found, 
we use a simple argument (see Section |3J) to prove the existence of the function. It 
is completely rigorous, and its application only uses Mathematica (or any similar 
package) to plot polynomials of degree up to « 50. Thus, the lower estimate for 
oil is computer-assisted in a more substantial way than the upper estimate (which 
is just "Mathematica-assisted" ) . 

In order to show that S2 < 2~ 1 / 2 — 4 x 10~ 6 and obtain neighborhoods of interior 
points in the connectedness loci, we use a covering argument inspired by |S| and 
|18j . At one point we need to check that a certain set is covered by the union of 
3 5 parallelograms, which we do using a computer. 

The paper is organized as follows. In Section 2 we provide the background on 
iterated function systems and discuss the relation between zeros of functions in B 
and connectedness of self-affine fractals. We then state our results. In Section 3 
we show how to find double roots close to a local minimum of a polynomial with 
certain properties. In Section 4 we establish the covering results and estimate «2- 
In Section 5 we prove the existence of gaps in Q2 and estimate 02- Section 6 is 
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devoted to some variants and generalizations. Section 7 contains proofs of several 
auxiliary results. 



An iterated function system (IFS) is a finite collection of (strict) contractions 
{fi, . . . , fm} on a complete metric space. Given such a system, there is a unique 
nonempty compact set E satisfying E = \Ji< m fi(E), called the attractor of the 
IFS, see [7J. We only consider IFS on M. of the form {TjX + bj}j< m where Tj are 
linear maps and bj £ W 1 . Their attractors are called self-affine. For the maps to 
be contractive (in some norm) it is necessary and sufficient that all the eigenvalues 
of Tj are less than 1 in absolute value. 

We investigate when attractors are connected in the simplest case m = 2. The 
following is well-known. 

Proposition 2.1 (see 0). The attractor E of an IFS {fi, $2} is connected if and 
only if fi(E) n f2(E) 7^ 0. (Of course, the "only if" direction is obvious.) 

For IFS of two affine maps there is a simple sufficient condition for connectedness. 
We can assume bi = without loss of generality, making a change of variable. 

Lemma 2.2 (folklore). Let {Tix, T2X + b} be an IFS of contracting affine maps, 
such that max{||Ti||, HT2H} < 1 in some operator norm, and | det(Ti)| + 1 det(T2)| > 
1. Then the attractor is connected. 

We include a proof for completeness (see Section 7). Next we specialize even 
more, assuming that T = T\ = T2, and state a criterion for connectedness in terms 
of zeros of power series. 

Let E = E(T,h) be the attractor of the IFS {Tx, Tx + b}, i.e., the unique 
nonempty compact set in W d satisfying 



2. Preliminaries on IFS and statement of results 



E = TEU(TE + b). 



(2.1) 



Observe that 



00 





n=0 



since the right-hand side is well-defined and satisfies (|2.1|) . 
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We can assume, without loss of generality, that b is a cyclic vector for T, that 
is, H := Span{T k h : k > 0} = M. d . Indeed, otherwise we can replace T by the 
restriction of T to H and consider the corresponding IFS on H. 

Combining Pr op osition 12 . 1 1 and (|2.2[) easily implies the following criterion, which 
is known, at least in special cases. Recall that B is defined in Q1.1JI . 

Proposition 2.3. Let T be a linear contraction with (possibly complex) eigenvalues 
Xj, f or 3 = l,---,m, having algebraic multipicities kj > 1, and geometric multi- 
plicities equal to one. Let b be a cyclic vector for T. Then E(T, b) is connected if 
and only if there exists f £ B such that 

f(X ) = ... = f^- l \X 3 ) = 0, j = l,...,m. (2.3) 

In particular, connectedness does not depend on b. 

Since this is a key statement relating connectnedness of self-affine sets to zeros 
of power series, we include a short proof in Section 7. Combining Proposition 12.31 
and Lemma 12.21 yields the following result. We denote by B the open unit disk in 
the complex plane. 

Corollary 2.4. Let A = {Ai, . . . , A m } C D and let kj = k(Xj) > 1 be such that for 
any nonreal A £ A, re have A £ A and k(X) = k(X). If YYJLi \^j\ kj — V^j then 
there exists f € B having zeros at Xj of multiplicity > k(Xj) for j = 1, . . . , m. 

In particular, we obtain that for any k > 1, every A G [2 _1 / fc ,l) is a zero of 
multiplicity > k for some power series in B. In |4, Section 3] it is asked whether 
there exist power series (or polynomials) with coefficients in {—1,0,1} having a 
k-th order root strictly inside the unit circle for arbitrary k. Corollary 12.41 answers 
the question for power series in a strong quantitative way, but the question for 
polynomials is much harder and remains open. 

From now on, we restrict ourselves to the case d = 2. Applying an invertible 
linear transformation as a conjugacy, we can assume without loss of generality that 
T is one of the following: 

where a, b, A, 7 are real, a 2 + b 2 < 1, b ^ 0, |A|, (7! < 1, and 7 / A. Each of the 
cases leads to a set which we call the connectedness locus for the corresponding 
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family of self-affine sets. Namely, we consider the sets 

M := {z = a + ib G B : 3/ G B, f(z) = 0}, 
M := {( 7 ,A)g(-1,1) 2 : 3/eB, /( 7 ) = /(A) = 0}, 
O := {AG (-1,1): 3/g£, /(A) = /'(A) = 0}. 

Thus, M,M,0 are essentially the sets of parameters for which the attractors in 
cases (i),(ii),(iii) respectively are connected. (It is natural to allow 7 = A in M and 
b = in A4 to ensure that the sets are relatively closed in D.) By Lemma 12.21 

M D M t := {A G D : |A| > 2~ 1/2 }, 
N D M:={(7,A)G (-1,1) 2 : | 7 A| > 1/2}, 
O D O t := {A G (-1,1) : |A| > 2' 1 / 2 }. 

We refer to A^i,A/i,Ot as "trivial parts" of the corresponding sets. 

The set A4 was studied by several authors, see [31 [21 03 El E EEU • In particular, 
Bousch [3] proved that M. is connected and locally connected. The set M has 
not been studied as much, although partial results are obtained in J7j, where it 
is shown that a large "chunk" of ftf is connected (all of M is conjectured to be 
connected). An approximation to AAn (0, l) 2 is depicted in Figure 1, with the non- 
trivial part shown in black. The picture is created with a program of C. Bandt; 
note that the visible disconnected pieces are a computing artefact. 




Figure 1 . Part of the connectedness locus M 
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By symmetry, we have O = O2 Pi (— O2). The set O2, defined in (|1.2j) . is our 
main object of study. Since B is a normal family, O2 is relatively closed in (0, 1), 
so there exists a% = minl^- 

Theorem 2.5. (i) a 2 G (0.6684755,0.6684757); 

(ii) disconnected. In fact, the intervals Ij, for j < 5, lie in distinct 

components of (0,1) \ n 2 where Ij = .668 + W~ 3 I'j, and I[ = (.478, .489), V 2 = 
(.632, .654), I' 3 = (1.282, 1.306), ^ = (1.327,1.333), I' 5 = (1.343, 1.352). 

This theorem is proved in Section 5, using results of Section 3 for part of the 
proof. 

Remark 2.6. Numerical evidence suggests that there are infinitely many compo- 
nents of (0, 1) \ ^2- We do not have a proof of that. The topological structure of 0,2 
appears to be very complicated. It can be proved rigorously that the five "gaps" 
above are the largest. 

Remark 2.7. Let 5 2 := sup((0, 1) \ 2 ). We have [2" 1 / 2 , 1) C 2 by Lemma l2~2l 
hence 0.2 < 2 -1 / 2 . We are able to show in fact that a.2 < 2 -1 / 2 — 4 x 10 -6 (see 
Theorem 12 . 1 01 below) . Our rigorous numerical results also yield 012 > .669355, and 
it seems that 012 < -67, but we do not have a proof. 

Returning to the sets M, J\f, we note that they are related to O. In fact, we 
have the following, denoting for Fct, Diag(F) := {(A, A) : A G F}. 

Lemma 2.8. (i) clos(.M \ R) n R C O; 

(ii) clos(AA \ Diag(R)) n Diag(R) C Diag(O). 

This follows by an easy compactness argument. See |151 Lemma 2.5] for the 
proof of (i); (ii) is proved in Section 7. 

As a consequence of Lemma 12.81 the pictures of Ai near the real axis reveal 
something about the structure of O. In fact, Figure 7 in served as an inspiration 
for our work. Lemma EH1 also explains an interesting feature of A4 and jV, namely 
the "antennas." The antenna for A4 (let us denote it by T(Ai)) is defined as 
the connected component of [|, 1) \ clos(7W \ R) containing ^ (there is obviously 
a symmetric antenna on the negative real axis); it was first discovered in [3] and 
studied in ^J^. Similarly, the antenna for M, denoted by T(M), may be defined 
as the connected component of Diag([^, 1)) \ clos(AT \ Diag(R)) containing (i, i). 
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Bandt noted that the "tip of the antenna" supr(.M) is the infimum of the set 
of double zeros of power series in B with infinitely many coefficients not equal to 
+1. It is not hard to show that supT(AA) is the infimum of the set of double zeros 
of power series in B with infinitely many coefficients not equal to —1. It seems very 
likely that supT(.M) = supT(AA) = but we do not know how to prove this. 
However, as a by-product of our investigation, we obtain the following corollary 
proved in Section 5. 

Corollary 2.9. We have sup r(.M), supT(AA) G (0.6684755,0.6684757). 

Our second main result concerns non-trivial interior points of the connectedness 
loci. In jHj and ^H] some chunks of interior points of M\Mt were found. Although 
numerical experimentation indicates that J\f\Nt and 0\Ot also have nonempty 
interior, this had not been proved rigorously before. 

Theorem 2.10. Let rj = 4 x 10" 6 . Then 



M D 


At = 


{(7,A)e(0,l) 2 


: | 7 _ 2 - 1 /2| ) | A _ 2 -i/2| <?? j ) 


(2.4) 


N D 


A 2 = 


|A G D : |A-2" 


- 1/2 | < 3^/4}, 


(2.5) 


O D 


A 3 = 


(A €(0,1): |A- 


-2-^1 < V \. 


(2.6) 



This theorem is proved in Section 5. The value of n is very small and clearly not 
optimal, but nevertheless it is an explicit constant. 

3. Existence of double roots I 

We denote by B n the subset of B consisting of polynomials having degree less 
than or equal to n. 

Let us say that (P, n, a, b) is good if P G B n , 0.5 < a < b < 1 (in reality, we will 
only consider 0.66 < a < b < 0.68), 

P(a) > a n+1 /(l -a), P(b) > b n+1 /(l - b), (3.1) 

P(x) > for all x G [a, b], and 

3 x G (a,b) : P(x) < x n+1 / (1 — x). (3.2) 

Lemma 3.1. Suppose that (P,n,a,b) is good. Let Q(x) = P{x) — x m where m is 
the minimal integer greater or equal to n + 1, such that Q(x) > on [a,b]. Then 
(Q, m, a, b) is good. 
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Proof. It is clear that Q £ B m . We easily check that 

Q(a) = P{a) -a m > a n /(l - a) - a m > a m+1 /(l - a), 

since a £ (0, 1), and similarly, Q{b) > 5 m+1 /(l — b). It remains to check the last 
condition. Either m = n + 1, in which case we note that 

Q(x) = P(x) - x n+1 < x n+1 /(l -x)- x n+1 = x n+2 /(l - x), 

or m > n + 1, in which case there exists t £ (a, b) such that P(t) — t m_1 < 0. Then 

Q(t) = P(t) -t m = (P(t) - t" 1 - 1 ) + {t m ~ l - t rn ) < t™- 1 - t m < f n+1 /(l - t), 

since t is greater than Clearly t ^ a,b, and the proof is complete. □ 

Corollary 3.2. Suppose that (P, n, a, b) is good. Then there exists f £ B such that 
P is the initial part of f , and f has a double zero in (a,b). 

Proof. Iterating Lemma 13. II we obtain a sequence of polynomials Qj £ B m . such 
that mj — > oo and Qj is the initial part of Qj+i for all j. Then Qj — > / G B 
uniformly on compact subsets of the unit disk in the complex plane. Since Qj are 
all positive on [a,b], we have that f(x) > on [a, b]. Since Qj(xj) < x™ 3 for some 
Xj S (a, b), we have that min^] / = 0. On the other hand, (|3.1|) implies that any 
function in B with initial part P is strictly positive at a and b. It follows that / 
has a zero in (a, b) of order at least two. □ 

Remark 3.3. The function f(x) = 1 + a j xJ obtained in Corollary 13 . 2l has the 

property that o, € {0, —1} for all j > n + 1 and there are infinitely many 0's and 
— l's. The only thing to check is that there are infinitely many 0's; the rest is obvious 
by construction in Lemma f3.ll Suppose that a, = — 1 for all j > N for some N £N. 
Then at some step in our construction we have good (Q,N,a,b), which implies 
Q(x) < x N+1 /{l - x) for some x £ (a, b). Then f(x) = Q{x) - x N+1 /(l - x) < 0, 
which is a contradiction. 

Corollary 13.21 is not very efficient numerically, since it allows us to find a double 
zero with an error of order 6 n / 2 . The next statement shows that this can be 
improved considerably. 

Corollary 3.4. Assume that (P,n,a,b) is good, b < .68, and n > 10. Further, 
suppose that P' has a zero at y £ (a, b), and 

P"{x)>C" for all x £ [a,b] (3.3) 
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Then there exists f G B such that P is the initial part of f and f has a double root 
in the interval (y — r],y + rj), where 

'■- 'Vy v ' <( ^ ( " +ir - 

Proof. From Corollary 13.21 we know that there exists / G B such that P is its 
initial part, and / has a double root in (a, b); let r be this root. Then P'(y) = 
f'(r) = whence, using the intermediate value theorem, 

\P'(r) - f'{r)\ = \P'{r) - P'{y)\ > C"\r - y\. (3.4) 

Note however that 

,p» - m < t *~ - ( - +i T;; )+r " +1 < 1+( ;r-l ( ; +v 

i=n+l ^ ' ^ ' 

(3.5) 

Combining (|3.4[) and (|3.5j) yields the corollary. □ 
Example 3.5. Consider 

P{x) = l - x 1 - x 2 - x 3 + x 4 + x 6 + x 7 + x 9 + x 10 + x 12 + x 13 

+ x 14 + x 16 + x 17 + x 18 + x 20 + x 21 + x 23 + ^24 + ^25 

+ x 26 + x 27 + x 28 + x 29 + x 31 + x 32 + x 33 + x 34 + x 35 

+ x 36 + x 37 _ x 38 + x 39 + ^40 + ^41 + x 42 + x 43 _ x 44 



+ x 45_ x 46_ x 47 + x 48 _ x 49_ x 50 GS5o> 

Let a = 0.668470, b = 0.668482. We claim that (P, 50, a, b) is good. We checked 
that P(x) > on [a, b] by plotting the graph of P on [a, b] in Mathematica; see 
Figure 
We have 

P(a)(l - a)a~ 51 » 1.03199, P(6)(l - b)b~ 51 » 1.06665, 

verifying (|3.1jl . On the other hand, 

p(r)(l - r)r~ 51 « 0.912958, for r = 0.6684756, 

verifying (|3.2j) . Thus, Corollary 13.21 applies, and we have a double zero in 
[0.668470,0.668482]. 
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Figure 2. Checking the existence of double roots: the polynomial 
P was obtained with the help of a C++ program. We plotted 
P(x)(l-x)/x 51 on the interval (0.668470,0.668482); it is clear from 
the picture that (P, 51, 0.668470, 0.668482) is good. 

Using Corollary 13.41 we obtain a more precise estimate. We have 

P(r + 2 x 10~ 8 ) - P(r) sn 8.89847 x 10~ 15 > 0, 

P(r - 2 x 10" 8 ) - P(r) » 2.04733 x 10~ 15 > 0, 

which implies that there exists y E (0.66847558,0.66847562) such that P'(y) = 0. 
We checked that P"(x) > 20 on [a,b] by plotting the graph of P" on [a, b] in 
Mathematica. Thus, Corollary l3~H applies. Since 506 50 • (4/20) < 2 x 10~ 8 , we 
obtain that there exists / E B with a double zero in (0.66847556,0.66847564). 

4. Existence of double roots II and connectedness loci 
Here we prove Theorem 12. 101 The proof is based on several lemmas. 

Lemma 4.1. Let T be a contracting linear map on IR 2 , and let b E M? be a vector 
such that T k h, k > 0, span M 2 . Denote the attractor of {Tx, Tx + b} by E. Let v 
be a point of the form 

k 

v = ^a i T- i b, 

i=l 
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where Oj G {—1,0, 1} and = 1. Assume that there exists a set U containing v 
and j G N such that 

Uc [j (T + «ib)o...o(r + «j-b)(tr). (4.1) 

ue{-i,o,i}3 

T/ten K is connected. 

Proof. There are similar results in e.g. [H], but we sketch a proof for completeness. 
Let E denote the attractor of the IFS {Tx— b, Tx, Tx+b}. Condition (|4.1|) implies 
that U is contained in E. In particular, v £ E. Recalling the form of v we get that 

k oo 

a-T^h = ^ c i Tih ' for some °i G {-!> °> !}• 

i=l i=0 

Applying T fe from the left on both sides we obtain a power series f £ B such that 
/(T)b = 0. Write f = f\ — fa, where /, has coefficients 0, 1 only. In particular, 
since at ^ 0, the power series /i and fi have different constant terms, hence 

A(T)b = / 2 (T)b G TE 1 n (TE + b). 

By Proposition 12.11 the set E is connected. □ 

Lemma 4.2. Let 

-( 2 r -(;)■ 

Lei afeo p = (—2.10, 0.20), q = (4.90,2.45). Denote by U the open parallelogram 
with the vertices ±p, ±q ; and denote by V the open parallelogram with the vertices 
±0.95p,±0.95q. Then 

v = T _1 b - T- 2 b - T~ 3 b - T~ 4 b + T~ 5 b ps (2.95837, 1.75736) G V, 

and 

Uc (J (T + uibJ-.-Cr + ttsbJCF). (4.2) 

ue{-i,o,i} 5 

Proof. This is the part of the proof that is computer assisted. The parallel- 
ograms U and V and the vector v were obtained through experimentation with 
Mathematica. 

The coordinates of v in the base {p, q} are (a,f3) (0.223,0.699). Since \a\ + 
\/3\ < 0.95, we get that v G V. We checked (|4,2j) rigorously using an algorithm 
which we now describe. 
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Figure 3. The parallelogram U is pictured in black. On the left, 
the gray figure is (TV — b) U (TV) U (TV + b); notice that U comes 
close to being covered by the iterates of V already in the first step. 
However, 5 iterations steps are needed in order to get a complete 
covering. This is illustrated by the picture on the right; there, the 
gray figure corresponds to the right-hand side in (|4.2|) . 

Given a parallelogram P, a collection of parallelograms Q and a depth n, we 
wrote a routine with the following pseudo-code: 

(i) If P is completely covered by some Q 6 Q, return true and exit (by 
convexity, it is enough to check whether all 4 vertices of P are contained 
in some Q). 

(ii) If n = return false and exit. 

(iii) Split P into 4 congruent pieces Pj. For each i, compute Qi, the family 
of parallelograms in Q intersecting Pj. 

(iv) If any Qi is empty, return false and exit. 

(v) For each i, run the routine with input Pj, Qi and n — 1. If all of them 
return true, then return true; otherwise, return false. 

Thus if the routine returns true then P is covered by the union of the family Q. 
This was indeed the case with the input P = U, n = 8 and 

Q= {(T + Ul h)---(T + u 5 b)(V) :ue {-1,0,1} 5 }, 

and this proves the lemma; Figure 01 depicts the situation graphically. □ 
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Proof of Theorem Let T, b, p, q, U, V be as in Lemma 14.21 Using Lemma 
14.21 we see that for sufficiently small 77, 

||r'-r||<?7 => 7c(T' + tj 1 b)o-..o(r' + u 5 b)(7). (4.3) 

The bulk of the proof will consist in obtaining an explicit value of 77 such that 1)4.3)) 
holds. 

Lemma 4.3. The condition \4-3j ) holds for n = 4 x 10 -6 . 

The proof of the lemma is straightforward, but technical, so we postpone it till 
Section 7. As mentioned above, the existence of such positive 77 is obvious. 
An immediate consequence of ()4.3jl and Lemmas 14. II and 14.21 is 

\\T — T\\ < rj ==> the attractor of {T'x, T'x + b} is connected. (4-4) 

It is convenient to use the £°° norm in M 2 and the operator matrix norm; recall 
that the latter is computed as the maximum of i 1 norms of the rows. Let 

Mi( 7 ,A) = h YJ' 7 /A, 
M 2 {r,e) = I _ £ J, p,e>0, 



M 3 (A) 



A 0.7 
A 



Note that Mi (7, A) is conjugate to the diagonal matrix with eigenvalues 7 and A; 
likewise, Ms(X) is conjugate to the standard Jordan block with eigenvalue A. Note 
also that 

||M 1 ( 7 ,A)-T|| = max(|2- 1 /2_ 7 | ) | 2 - 1 /2_A|); 
||M 3 (7,A)-r|| = p-^-AI. 

From this, Lemma EOl and ()4.4|) we obtain ()2.4|) and l|2.6jl . 

It remains to consider the complex eigenvalue case. The matrix M2(p,e) has 
eigenvalues p±iy/0.7e. Let A be a non-real complex number such that |A — 2 -1 / 2 | < 
(3/4)77; without loss of generality assume that Im(A) > 0, and write A = p + i^/oTfe. 
We have that \p - T" l l 2 \ < (3/4)77 and 

V 7 ^ < (3/4)77 => £< 1R ^ n 7 77 2 < (1/4)77. 

lb x U.7 
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Therefore \p — 2~ 1//2 | + e < r], and from this we conclude that ||M2(p, e) — T\\ < r/. 
Invoking (|4.4|) again, this shows that (|2.5|) is verified, which completes the proof. □ 

Remark 4.4. Although the proof of Theorem l2.1(Jl is inspired by analogous results 
for self-similar sets which appeared in [S] and JT5J, the more complicated geometry 
of self-affine sets introduce some additional difficulties. For example, the use of a 
computer to first find the sets U, V and the vector v, and then for checking that the 
covering property holds, becomes essential (in the self-similar case, some attractors 
are actually rectangles, which allows to do a purely algebraic analysis in certain 
region, see |18j). The chunks of interior points of A4 \ Mt which were found in 
those papers are away from the real line. 

Remark 4.5. Let i^ 7i A be the attractor of {T 7 a, T 7 a + (1, 1)}, where T^ : \ is a diag- 
onal map with eigenvalues 7, A. In |13j the almost sure Hausdorff dimension of K^ x 
was found in the region (0, l) 2 \A/f. The result was new only in the region M\Mt', 
hence Theorem 12. 101 makes that result effective, by showing that C2(N\Nt) > 0. 

Remark 4.6. The value of r/ found in Lemma f4.3l is extremely small, but graphical 
experimentation suggests that the same covering argument, even with the same 
covering, works for a large range of parameters. We believe that it should be 
possible to extend the result to show that actually (0.7, 2 -1 / 2 ) C (we recall that 
computer results suggest that indeed (0.67, 2 -1 / 2 ) C CI2, but this appears harder 
to prove rigorously). 



5. Existence of gaps 

We describe an algorithm which we use to rigorously prove the existence of gaps 
in the set O2. As a by-product, the results of Section El yield bounds on how large 
the gaps may be. In practice, we are able to obtain an accurate description of the 
set Q2 up to an error of 10 -8 . 

Our algorithm is based on Bandt's algorithm to study the Mandelbrot set for 
pairs of linear maps [Q. The idea is the following: assume that / S B has a double 
zero in (a, b) C (1/2, 1), and let P n G B n be the initial part of / up to exponent n. 
Then, letting r be the double root, 
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00 un+1 

\P n (r)\ = \P n (r)-f(r)\< £ r l < - 



i=n+l 



V 



\K(r)\ = |P»-/V)I< E »'-< (1+ '' < !" 2 '" r . 

i=n+l ^ ' 

Hence, using that a < r < b and the intermediate value theorem, 

\P n (b)\ < \P n (r)\ + ||f£b~(r,6)(&-r) < ^-^ + *~ ° ; 

1^(6)1 < 1^(01 + ||^||xo.(r,6?(& - r) < (1 + ( i ( ^ )a * ))y> + 

Conditions (|5.1|) are easily checkable. If for some n and (a, 6), at least one of 
them fails for all P € B n , then (a, 6) fl Q2 = 0- 

We discuss some of the details of our implementation. Full C++ code is available 
upon request. 

We used floating point arithmetic in order to get practical performance (in prin- 
ciple, one could use exact rational arithmetic instead). To keep the algorithm 
rigorous we computed the theoretical floating-point error of our calculations and 
added a corresponding error term to the inequalities (|5.1|) . In order to reduce the 
number of arithmetic operations and avoid very small numbers and the consequent 
loss of precision we multiplied both sides of ()5.1|) by b~ n . Thus the numerical checks 
that we used are 

L-nnwLM b (b-a)b- n 

b- n \P'(b)\ < 1 + n ^- 6 ) 1 (P~a)b- 1 „ (52 ) 
b l^nWI < (1 _ 6)2 + 2(1 _ 6) 3 +T I- ^ 

Here 77 is the error term. In practice taking r\ = 10~ 14 suffices. This was calcu- 
lated based on the IEEE floating-point standard. We used the MinGW compiler 
on a Windows XP platform. 

In order to make the algorithm efficient we exploited the tree structure of B n . 
The basic routine takes as arguments an interval (a, 6), a depth d, and a polyomial 
P £ B n . The routine returns a boolean value, which we denote by C((a,b),d, P). 
This value indicates whether the inequalities (|5.2|) are verified for at least one 
polynomial of degree at most n + d with initial part P. Hence if C((a, b),d, 1) 
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returns false for some d we must have (a, b) n fl2 = 0- The structure of the routine 
is as follows: 

(i) Check (|5.2[) for P and the interval (a,b). If any of the inequalities fails 
to hold, return false and exit. 

(ii) If d = 0, then return true and exit. 

(hi) For i = —1,0 and 1, run C((a,b),d — l,P(x) + ix n+1 ). If any of these 

returns true, then return true and exit, 
(iv) Return false. 

Assuming that C((a, 6), d, 1) returns true, this routine easily produces a poly- 
nomial of degree <i for which ()5.2|) holds (this is done by keeping track of which i 
produces a true in Step [5] of the routine). This polynomial can in turn be used to 
show that there is a double root near b using the results of Section |21 

Note that in order to use a large depth one needs to run the algorithm on very 
small intervals, due to the presence of the term (6 — a)b~ n in the right-hand side of 
(|5.2j) . Even on a standard desktop PC, it took less than 3 hours to scan the interval 
(0.66847, 0.66936) for gaps using a grid of 10 7 subintervals and running the main 
procedure on each. Figure summarizes our findings. We plotted the logarithms 
of the lengths of the 60 gaps we found, as well as the lengths of the complementary 
intervals. Although we do not know how to prove it, we believe that at least the 
few largest of those pieces which the algorithm did not rule out contain intervals 
of the set O2 of approximately the same length. 

We finish this section explaining how we verified the assertions of Theorem 12.51 

Proof of Theorem \2. 51 . We used the algorithm to show that (0.5,0.6684754) and 
Ij (1 < j < 5) are contained in (0, l)\f22- We did not need to go deeper than 50 
iterations for this. For example, we subdivided the interval I\ into 110 intervals 
(ai,bi) of length 10 -6 and ran C((aj, 6j), 40, 1) for each i. The routine returned a 
false value for all, which implies that I± £ (0, 1)\^2- 

It follows from Example 13.51 that 02 < 0.6684757. To check that the intervals 
Ij lie in different connected components of (0, 1)\^2 we used the algorithm to 
produce suitable polynomials, and then followed the scheme of Example 13.51 to 
show the existence of points 0,(1 < j < 4), 9j £ (max Ij, min Ij+i) ^2- 

The following table summarizes our findings. The symbols on the right represent 
the coefficients of the polynomials Pj\ corresponds to the coefficient —1. The 
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length of n-th gap 



8 - 




i i i 

i 

i 

■ ■ 

6 - 



10 20 30 40 50 

Figure 4. This figure shows the distribution of the largest gaps in 
the set r^2- The negative of the base 10 logarithms of the length of 
the gaps are contained in the vertical segments (this is a rigorous 
consequence of the algorithm applied) ; heuristically, we believe that 
the actual lengths are closer to the upper end of the intervals. 



numbers 9j are within 10 7 of a double root; more precisely, for each j there exist 
fj £ B with initial part Pj, and 0j a double root of fj such that \9j — 9'A < 10~ 7 . □ 



j 






1 


0.668550 


10001O11O11O111O111110O1111100100O11111111111100100 


2 


0.668900 


l000lOllOllOlllll0lll00ll00ll0lll0lO00O0O0000O0O000 


3 


0.669310 


I000IOIIOIIIOOIIIOIOIIIIOIIIIII0OOO00000I0000IOOOII 


4 


0.669336 


I000IOIIOIIIOOIIIOIIOOIIIIIII0OIIOIOIIII0IIIIIII00O 



Proof of Corollary \2.fA The double zero y obtained in Example 13.51 is a root 
of / E B which has infinitely many 0's among the coefficients, by Remark 13.31 
Replacing a,j = by +1 for j large yields a function in B with a complex zero close 
to y. On the other hand, replacing a,j = by — 1 for j large yields a function in B 
with two real zeros close to y. Thus, y £ clos(.M \R) and (y, y) £ clos(AA\Diag(M)). 
In view of Lemma 12.81 this implies both "tips of antennas" belong to (ct2,y), a and 
the claim follows from Theorem 12.51 □ 
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6. Variants and generalizations 

The algorithm described in the previous section adapts without difficulty to more 
general settings; we consider some examples below. 

6.1. The set of triple roots. Here we restrict ourselves to the family B, but 
consider higher order roots. Geometrically, the set of roots of multiplicity n corre- 
sponds to the connectedness locus of self-affine sets associated to Jordan blocks of 
order n; recall Proposition 12.31 We denote this set by Q n . 

The computer algorithm extends in a straightforward way to higher multiplicity 
roots. Indeed, it is enough to replace (|5,1|) by the set of tests 

\P i Hb)\<H%\b) + ( , + *- a _ b)i+2 , 0<*<n, 

where H n (x) = x n /(l — x). The algorithm yields intervals in (0, l)\f2 n ; in partic- 
ular, it gives lower bounds on a n = min$7 n . We remark, however, that in practice 
the program becomes very slow for n > 4. For n = 3, we have the following result. 

Proposition 6.1. (i) a 3 = minf2 3 > 0.743, 

(ii) (0.746,0.7465) C (0,l)\n 3 , 

(iii) (2- 1 /3 ) i )c f] 3 . 

Remark on proof: (i) and (ii) are direct applications of the algorithm, while (iii) 
follows from Lemma l2~2l (Note that 2~ 1 / 3 « 0.7937). □ 

Remark 6.2. Numerical experimentation suggests that a 3 G (0.743,0.744) and, 
in particular, f2 3 is disconnected. However, the techniques of section |3] do not seem 
to apply, so this remains a conjecture. 

It is interesting to compare our results with those of 0], where multiple roots of 
the following family were considered: B = {l + a i x% : a* S [ — 1,1]} . Let [3 n 
be the smallest root of multiplicity (at least) n of some / E B. The values of (3 n 
for n < 27 were computed in 0]; in particular, ~ 0.64914 and /? 3 w 0.72788. 
Observe that 02— 02 > 0.01934 and a 3 — /? 3 > 0.0151. Thus going from a continuous 
to a discrete set of coefficients does have a substantial impact in the set of multiple 
roots. 



6.2. The set of double zeros with coefficients 0,±1,±2. We can generalize 
the set O2 in another direction by enlarging the set of allowed coefficients. For 
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concreteness, we will work with the coefficient set {0, ±1, ±2}; specifically, let 

oo 

B' = {l + a ^ ■■ a i G {°> -!> li -2= 2 }}- 

i=i 

Denote by Q' 2 the set of double zeros in (0, 1) of elements of B' . It turns out that 
a specific power series plays a very special role in the study of this set. Let 

Q(x) = 1 - 2x - 2x 2 + > 2x l = 1 - 2x - 2x 2 + . 

l-x 

Note that 1/2 is a double root of Q(x). In fact, more is true: Q(x) is a so-called 
(*)-function for the class B' on the interval (0, xo) for all xq < 1/2; see ^H] for the 
relevant definitions and proofs. A consequence of this is that minf^ is precisely 
1/2. 

Here we prove that 1/2 is actually an isolated point of f^; this is due to the 
special form of the function Q(x), and in particular the fact that all but finitely 
many coefficients are +2. More precisely, we have the following result: 

Proposition 6.3. min(fi£\{l/2}) G (0.5436,0.5438). 

Before proving the proposition we remark that the set-up of Section |21 works here 
with minor modifications. In this context, we say that (P, n, a, b) is good if P E B' n 
(the family of polynomials in B' of degree at most n), 0.5 < a < b < 1, 

P(a) > 2a n+1 /(l - a), P(b) > 2b n+1 /(l - b), 

P(x) > for all x G [a, b], and 

3xe{a,b): P(x) < 2x n+1 /(l - x). 

The proofs of Section |21 apply almost verbatim. In particular, if (P, n, a, b) is good 
then there exists a sequence %\ < %i < . . . such that f(x) = P(x) — YljfLi < ^ x%3 ^ as 
a double root in (a, b). 

Proof of Provosition \6.'J\ We will use the following result, which follows from a 
modification of the proof of Theorem 2 in [I] (see |13j for a complete proof): if / 
is a power series with coefficients in [—1,1], and ai, . . . are complex roots of / 
in the unit disk, counted with multiplicity, then 

/ 1 \ - fc / 2 

K ... a k \ > I 1 + -J (fc+l)- 1 / 2 . (6.1) 
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We will use this result with k = 3. A standard application of the algorithm shows 
that / 6 B' may have roots in (0.5,0.51) only if it starts with 1 — 2x — 2x 2 . Then 
it follows from the definition of Q that f{x) < Q(x) for all x > 0. Since /(0) = 1 
and /(1/2) < Q(l/2) = 0, / has a root in the interval (0, 1/2). Suppose that a is 
a double root of / in (1/2, 1). We obtain from ()6.1|) that 



the interval (0.5,0.51), whence (0.5,0.51) C (0, 1)\0^. 

A standard application of the algorithm shows also that (0.51, 0.5436) C (0, 1)\02- 
Finally, let P be the polynomial of degree 26 with coefficients 

(1,-2,-1,1,1,1,2,1,1,2,1,1,2,1,1,2,1,2,1,1,2,-2,-2,-1,-2,2,-1) 

We checked with Mathematica that (P, 27, 0.5436, 0.5438) is good. This implies 
that minf^ £ (0.5436,0.5438), completing the proof. □ 

Remark 6.4. The set Q' 2 \{l/2} seems to be connected (i.e. no "gaps" appear 
when running the program), but we do not have a proof of this. Still, if what 
the numerical experimentation suggests holds true, then the sets ^2 and £1' 2 have 
strikingly different topological structure. 



Proof of Lemma \2."A Let E be the attractor of the IFS, that is, E = T\E + 
(T2E + b). Let || • || be the norm in W d such that ||2~ix|| < r[|x[|, i = 1,2, for some 
r G (0, 1) and all x G W*. Denote by F e the e-neighborhood of a set F C M. d in this 
norm. Then (TiF) £ D TiF £ / r . 

Suppose that E is disconnected. Then T\E n (T2E + b) = by Proposition 12.11 
and since these are compact sets, we can find e > such that {T\E) £ P\ (T2-E' + b) e = 
0. Then E £ = {T\E) £ U (T2E + b) e is a disjoint union, so 




whence a > (6-V/3/16) 1 / 2 > 0.8. We conclude that / cannot have double roots in 



7. Remaining proofs 



= ^((TrE^ + C^E),) 
> C d {T x Ez/ r ) + C d {T 2 E £/r ) 
= (| det(Ti)| + I det(T 2 )|) • C d (E e/r ) > C d (E e/r ). 



This is a contradiction, since E £ / r \ E £ has positive Lebesgue measure. 



□ 
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Proof of Proposition^^ If E = E(T, b) is connected, then TE n (TE + b) ^ 0. 
In view of (|2.2|) . we obtain that there exist {0, 1} sequences {a n }o° an d { a n}o° sucn 
that a = l,ag = 0, and X^L ( a n - a' n )T n b = 0. Denoting f(x) = Y^=o( a n ~ 
a' n )x n we get f £ B and f(T)h = 0. Now let us write b = Y^jLi c j e j where 
ej G Ker(T — Xjl) k i. Since b is a cyclic vector for T, we have Cj ^ and £" 
Ker(T - Xjl)^' 1 for j < m. Then /(T)b = implies that f(T)ej = for all 
j < m. We have 

/(T) ej = /(A,)e, + /'(A,)(T - A,!)^- + . . . + f^~ l \X 3 )(T - A ; /) ; ' 

Since the vectors {(T — Xjl) e ej : I = 1, . . . , kj — 1} are linearly independent, ()2.3|) 
follows. 

Conversely, if /(x) = 1 + E^°=i h nX n £ £ satisfies (JUSJ), then /(T)b = for all b. 
Writing b n = a n — a' n for some a n , a' n £ {0, 1}, we obtain that TE n (TE + b) ^ 0, 
hence £7 is connected by Proposition 12.11 □ 

Proof of Lemma W!\ (ii) Suppose that A £ (—1, 1) and (A, A) is such that there 
exists a sequence (7™, A n ) £ M, with 7„ < A n , converging to (A, A). Then there are 
power series f n £ N such that / n (A n ) = f{j n ) = 0- By compactness, passing to 
a subsequence, we can assume that f n —>f£B coefficientwise. Then fn k \x) — ► 
/( fc )(A) for all fe > 0. By Lemma 1X21 |A| < 2~ 1 /2. Let C fc = max{|/( fc )(x)| : / G 
,8, |x| < 2" 1 / 2 }, which is finite (and easy to compute explicitly). Then we have for 
n sufficiently large: 

|/n(A)| = |/n(A) " fn{X n )\ < C X \X - X n \ - 0, n - OO. 

Next, there exists i n G (7„, A n ) such that f n {t n ) = and we have for n sufficiently 
large: 

l/n(A)| = |/4(A) - < C 2 |A - i n | - 0, 

It follows that f(X) = /'(A) = 0, as desired. □ 

Proof of Lemma In the following calculations ry will be a fixed positive 

number, to be determined later. We use the £°° norm on M 2 . 
Given a word u G {—1,0, l} 5 let 



T u = (T + mb)o...o(T + u 5 b), 
T' u = (r' + «ib)o---o(T' + 'U5b). 
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Observe that 



8=1 

(aft -1 * = (rO- 6 x-(r')- 5 ^« i (ro < - 1 bj. 

Hence 

5 5 

CO _1 r«x = ((T')- 5 T 5 ) x + (T')~ 5 ^ u^b - (T')- 5 ^ Ui(T'y~ 

i=l i=l 

=: 5x + d, 

where 

5 

5 = (T')- r °T 5 , d = (T')~ 5 J^T* -1 - (T') i_1 )b. 
Observe that (for r? < 10" 



i=l 

1-2 A 



H^-rfe _ (r') fc || < ^ [|T J '(T') k ~i — 7 1 i~ 1 (7 1 ') fc ~i+ 1 | 
i=l 

= J2l|T 3 '- 1 (r-r')(T0 fc ~ i ~ 1 || 

i=i 

< fcrymaxdlTlUlT'll)^ 1 < fc (1.5)*- 1 r?. 
Hence if it! = (T') 5 — T 5 then ||i?|| < 5 (1-5) 4 r\ < 26rj, and we can estimate 

H(T')- 5 || = || (t 5 (i -t^r))' 1 || < ||(/-r- B i?)- 1 ||||r- 5 || 



< r- 5 ||^(||r- 5 |||| J R||)^ = r l|r 



- ||T-B||[|/E|| 
33.6 

< < 34, 

1 - 33.6 x 26 7/ 

as long as rj < 10~ 5 . Therefore, 

||5 - J|| = IKT')- 5 ^!! < ll(?T 5 ||||i?|| < 34 x 267? = 884 r,. 
We have by l|72|) . 

4 

||d|| < IKT') -5 !) ^ j (1.5) J ' _1 ?7||b|| < 34 x 24.25 X] < 825 77. 
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Note that 

max{||x|| : x € V} = 0.95max(||p||, ||q||) = 0.95 x 4.9 < 4.7. 

It follows from the previous estimates that for x £ V, 

\\Sx - x|| < \\S - J||||x|| < 884 x 4.7?? < 4155??. 

In particular, this implies that SV C Vg, where 5 = 4155 rj, and V$ denotes the 5 
neighborhood of V. Recalling (|7.1|) we get that 

(Ti)- l T u V C {V s ) m = V 5+m C Vp, (7.5) 

where 8' = 5 X 10 3 rj, since ||d|| < 825 r/ by (|7.4j) . 

Let W(r) = {(x,y) : \x\ + \y\ < r} and let M be the matrix with the columns 
p, q. Note that U = MW(1) and V = MW(0.95). Note also that 

dist(iy(0.95),M 2 \iy(l)) = 0.05/2 = 2.5 x 10~ 2 , 

where dist(-, •) denotes the distance induced by the £°° norm. An easy calculation 
yields ||M _1 || = 1.2. Therefore, 

, ,. r ^rr, dist(W(0.95),R 2 \W(l)) 2 .„ , 

dist(F,M 2 \lT) > i — v HM- 1 !! — > X ^ ^ 

It follows from ()7.5|) and (|7.6|) that if r\ is so small that 5' < 2 x 10~ 2 , then 

(r^)- 1 ^^ C U. (7.7) 

Since 5' = 5 X 10 3 r/, this will be the case for rj = 4 x 10~ 6 . From now on, we will 
fix this value of r\ (since rj < 10~ 5 , the previous calculations apply). 

Let A be the "multiplication by 0.95" map, so that V = AU. By applying A to 
both sides of (|4.2j) we get 

V C A j |J r tt (vo I = |J T U AV. 
\ug{-i,o,i} 5 / «e{-i,o,i} 5 

However, we deduce from (|77J) that T U V C T^A^V, or T' U V D T U AV . Combining 
this with the last displayed formula, we conclude that (|4,3[) holds. □ 
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